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PERIODICITY  AND  PECULIARITY  IN  120 
FIRST  AND  SECOND  ROW  DIATOMIC  MOLECULES 


Alexander  I.  Boldyrev,  Nick  Gonzales  and  Jack  Simons 
Department  of  Chemistry,  The  University  of  Utah,  Salt  Lake  City,  Utah  84112 

Abstract 

The  ground  and  very  low-lying  excited  states  of  all  120  first  and  second 
row  diatomic  molecules  are  surveyed.  Three  quarters  of  these  molecules  have 
had  their  ground  state  term  symbols  reliably  experimentally  determined. 
However,  one  quarter  remain  predicted  only  theoretically.  For  all  120  species, 
the  best  available  experimental  (where  known)  and  theoretical  values  for  the 
dissociation  energies  to  ground-state  atoms  are  also  presented.  The  Aufbau 
principle,  combined  with  standard  energy  ordering  for  the  valence  molecular 
orbitals,  is  able  to  properly  account  for  the  ground-state  term  symbols  of  all  but 
twenty  of  the  diatomics  studies.  The  twenty  exceptions  produce  higher  than 
expected  ground  state  spin  multiplicity  and  arise  when  there  are  4-5  or  7-8 
valence  electrons  and  group  3,  4,  or  5  (but  not  group  6  or  7)  atoms  are  involved. 


I.  INTRODUCTION 

One  might  expect  that  essentially  all  of  the  15x16/2=120  diatomic 
molecules  comprised  of  first  (H,  Li,. . .  F)  and  second  (Na, . .  .Cl)  row  atoms  have 
been  thoroughly  studied  to  the  extent  that  their  ground  electronic  states  and 
corresponding  bond  lengths  (Re)  and  dissociation  energies  (De)  are  well 
established.  However,  such  is  not  the  case;  in  Fig.  1  those  diatomics  for  which 
even  the  ground  electronic  states  have  not  been  so  characterized  are  displayed 
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in  burgundy.  In  purple  are  shown  the  diatomics  whose  ground  electronic  states 
are  reasonably  well  characterized.  Most  of  the  experimental  data  used  to  create 
Fig.  1  was  taken  from  the  monograph  of  Huber  and  Herzberg,''  although  several 
species'  properties  were  obtained  from  more  recent  sources.2-9  It  probably 
surprises  most  students  of  chemistry  to  learn  that  more  than  one  quarter  of  all 
the  diatomic  molecules  formed  by  combining  pairs  of  first  or  second  row  atoms 
have  yet  to  be  experimentally  characterized.  Many  of  the  uncharacterized 
diatomic  molecules  are  very  reactive  intermediates  with  unpaired  electrons  or 
unsaturated  valences  of  one  or  both  atoms,  which  therefore  can  exist  and  be 
studied  only  under  special  conditions. 

(Figure  1  near  here  please) 


In  this  article,  we  consider  the  electronic  structures  of  the  ground  and 
low-lying  excited  states  of  diatomic  molecules  composed  of  atoms  from  the  first 
and  second  rows,  including  the  corresponding  hydrides  but  excluding  rare-gas 
containing  species.  We  emphasize  (i)  species  that  have  yet  to  be  studied 
experimentally,  (ii)  species  whose  ground  states  do  not  involve  maximal  double 
orbital  occupancy,  (iii)  trends  and  exceptions  to  trends  in  the  spin  multiplicity  of 
ground  states. 

Sophisticated  ab  initio  techniques  were  applied  to  many  of  the  thirty- 
three  experimentally  uncharacterized  diatoms  shown  in  Fig.  1  burgundy.  In 
particular,  the  following  twenty-three  have  been  studied  in  earlier  theoretical 
works:  LiB.io  LiC.^Oc.ii  LiN.i0c-i2  LiMg,i3  kiAl.iOc  LiSi,iOc-ii  Lip,i4  LiS,‘'5 
BeN,'>6  Be^i^  BeNa,''3  BeA^s  BeSi,''^^  BeP,^'*  BP,i^  NaMg,i3  NaSi,iOc 
NaP,""^  NaS,‘'5  MgSi,''®^  MgP.i^  AIP,''^  and  SiP.''^  In  the  present  work,  we 
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present  our  new  results  on  the  remaining  ten  diatomic  molecules:  BeB,  NaB, 
NaC,  NaN,  MgB,  MgC,  MgN,  AIB,  NaAt  and  MgAI  and  we  repeat  high  level 
calculations  on  several  of  the  other  twenty-three  molecules  for  which  the  ground 
state  has  not  yet  been  identified  with  certainty.  In  addition,  we  attempt  to 
examine  patterns  in  ground-state  spin  multiplicity  for  the  species  in  Fig.  1,  in 
particular  noting  circumstances  where  ground  states  with  higher  than  expected 
spin  multiplicities  occur. 

II.  COMPUTATIONAL  DETAILS 


The  bond  lengths  and  harmonic  vibrational  frequencies  of  the  ten 
diatomics  (BeB,  NaB,  NaC,  NaN,  MgB,  MgC,  MgN,  AIB,  NaAl  and  MgAI)  for 
which  new  data  is  presented  here  were  optimized  using  analytical  gradients'*® 
and  polarized  split-valence  basis  sets  of  6-31 1+G*^®  quality  at  the  correlated 
MP2(full)  level  (UMP2(full)  for  open-shell  systems)  of  theory.  The  resulting  MP2 
(full)/6-311+G*  equilibrium  geometries  were  then  used  to  further  evaluate 
electron  correlation  corrections,  in  the  frozen-core  approximation,  by  Moller- 
Plesset  perturbation  theory  to  full  fourth  order20  and  by  the  (U)QCISD(T) 
method2''  using  the  6-31 1-»-G(2df)  basis  sets  for  Li  to  F  and  Na  to  Cl.  The  UHF 
wave  functions  for  open-shell  systems  were  spin-projected  to  produce  pure 
spectroscopic  states  (PUHF,  PMP2,  PMP3  and  PMP4).22  The  geometries  of 
selected  low-lying  excited  electronic  states  were  optimized  then  at  the 
QCISD(T)/6-311+G(2df)  level.  All  calculations  were  carried  out  with  the 
GAUSSIAN  9223  suite  of  programs  unless  otherwise  specified,  and  core 
orbitals  were  kept  frozen  in  all  correlated  calculations. 

The  geometries  and  harmonic  vibrational  frequencies  of  several 
diatomics  that  have  two  electronic  states  within  5  kcal/mol  of  one  another  were 
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also  studied  using  the  complete  active  space  self-consistent  field  -  multi¬ 
reference  configuration  interaction  method,  including  all  single  and  double 
excitations  from  the  reference  configurations  (CASSCF-MRClSD)24-26  in 
diatomic  molecules  treated  by  CASSCF  and  CASSCF-MRCISD  calculations  in 
C2v  symmetry,  the  first  (ai)  representation  contains  both  o  and  5  orbitals,  the 
second  (bi)  and  third  (b2)  contain  jtx  and  Tty  orbitals,  respectively,  and  the  forth 
(a2)  contains  5-orbitals.  Because  we  used  different  active  spaces  for  different 
molecules,  we  present  the  details  of  each  calculation  when  discussed 
specifically  later.  For  the  CASSCF-MRCISD  calculations,  we  used  the  very 
large  ANO  basis  set  of  Widmark  et.al.27.  These  calculations  were  performed 
using  the  MOLCAS-2  program. 28 

The  MgAI  and  BAl  molecules  were  also  studied  at  the  MCSCF  level 
using  Dunnings'  augmented  correlation  consistent  polarized  valence  basis  sets 
(aug-cc-pVDZ)29  and  the  GAMESS^o  program.  Details  of  the  configuration 
spaces  employed  will  be  given  in  the  description  of  the  properties  of  these 
molecules. 

The  ground  and  low-lying  state  results  for  the  ten  newly  examined 
diatoms  as  well  as  for  five  others  examined  by  earlier  workers  and  reexamined 
here  are  summarized  in  Tables  l-XV. 

III.  RESULTS  AND  DISCUSSIONS. 

A.  The  Ten  As  Yet  Uncharacterized  Diatomics 

BeB.  Assuming  doubly  occupancy  for  the  la  valence  orbital,  the  three 
other  valence  electrons  may  be  distributed  throughout  the  lowest  2a,  3a  and  l  n 
valence  orbitals  giving  five  possible  occupancies;  1a22a2l7c^  1a22a23a^ 
1a22aM;t^3a^  1a22aM7r2  and  1a2l7r3  which  lead  to  several  low-lying 
electronic  states:  ^Ur  (1a22a2ln:i),  (1a22a23a‘'),  ^Tlr  and  ^Or 
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(1o22aMrt'*3o^).  and  {^a^2a'^^K^)  and  Preliminary 

calculations  at  MP2(full)/6-31 1+G*  level  were  carried  out  for  ^Ur,  "^Flr.  ^X‘, 
and  ^rii  states  with  the  ^Flr  state  proving  to  be  the  lowest  (Table  I).  However, 
when  larger  basis  sets  and 

(Table  I  near  here  please) 

more  sophisticated  correlation  methods  (QCISD(T)/6-311+G(2df))  were  used, 
the  ^Ur  state  was  predicted  to  be  the  ground  state,  and  the  2X+  state  to  be  the 
first  excited  state.  The  ^Or  and  states  are  the  next  excited  states.  For  BeB 
we  are  confident  that  the  state  is  the  ground  electronic  state  because  (i)  all 
four  low-lying  states  have  small  spin-contaminations,  (ii)  the  energy  difference 
between  the  ground  electronic  state  and  the  first  excited  state  is  8-9  kcal/mol, 
and  (iii)  the  relative  energy  difference  between  PMP4  and  QCISD(T)  is  only  1.3 
kcal/mol,  less  than  the  first  excitation  energy.  The  calculated  dissociation 
energy  (De)  of  BeB  (2nr)  is  1.57  eV  at  the  QCISD(T)/6-31 1-hG(2df)  level. 

NaB.  Given  a  total  of  four  valence  electrons  and  assuming  double 
occupancy  for  the  la  orbital,  the  only  three  low-energy  configurations  involve 
1a22a2,  Ia22a^l7:^  and  Io^Ik^  occupancies  all  of  which  have  been  studied. 
The  results  appear  in  Table  II  were  we  find  the  high  spin  state 

(Table  II  near  here  please) 

to  be  the  ground  state  and  the  low  spin  (Ia22a2)  state  to  be  the  lowest  lying 
excited  state  at  both  the  PMP4  and  QCISD(T)  levels. 
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Because  the  adiabatic  X^FTr  — » excitation  energy  is  rather  small,  3.1 
kcal/mol,  we  carried  out  calculations  using  the  CASSCF-MRCISD(Q)  method 
and  the  (l7s12p5d4f/7s5p3d2f)Na  +  (l4s9p4d3d/7s5p3d2f)B  basis  set.  The 
CASSCF  expansion  included  152  and  160  (^Flr )  configurations,  which 
represent  all  possible  excitations  of  all  four  valence  electrons  among  (4,2,2, 1) 
active  orbitals,  respectively,  of  ai,  bi,  b2  and  a2  representations  of  C2v 
symmetry.  All  single-  and  double-excitations  from  these  152  and  160  CASSCF 
configurations  were  then  included  in  the  MRCISD  calculations  giving  the 
86,114  CS-*-)  and  117,475  (^flr )  MRCISD  configurations.  Davidson's 
corrections  to  the  MRCISD  energies  were  very  small  and  as  a  result  the 
CASSCF-MRCISD  and  CASSCF-MRCISD(Q)  level  are  very  close  to  each  other 
(Table  I!).  Moreover,  the  ^Flr  state  is  the  ground  state  and  is  the  first  excited 
state  at  both  the  CASSCF-MRCISD  and  CASSCF-MRCISD(Q)  levels.  The 
ground  electronic  state  is  well  represented  by  the  single  Hartree-Fock  electronic 
configuration  (CHF(1<y22a‘' 1rt‘')=0.9525)  which  is  the  only  configuration  in  the 
MRCISD  expansion  that  has  a  coefficient  larger  than  0.15.  The  occupancies  of 
the  natural  MOs  for  the  ^Or  state  are: 

1  892a0.983<j0.03i  ;t^0.97i  ny0.062KxO  022Ky0  01 . 

The  ^1+  state  is  also  well  represented  by  the  Hartree-Fock  wave  function 
(Chf(1o22o^)=0.9034).  However,  in  this  case,  one  non-Hartree-Fock 
configuration  has  an  expansion  coefficient  larger  than  0.15:  C(1o22o®3o2)=- 
0.2707.  The  occupancies  of  the  natural  MOs  for  the  state  are: 

I  ^1.882<y1 -743^^0.171  I  ;Cy0.07  271x0  0^27tyO  O^ 

We  are  confident  that  ^Tlr  is  the  ground  electronic  state  for  NaB  because 
we  have  very  good  agreement  between  the  adiabatic  excitation 

energies  with  all  four  sophisticated  ab  initio  methods.  The  valence  isoelectronic 
LiB  molecule  also  has  a  XSFlr  ground  electronic  state,  however  other  valence 


isoelectronic  molecules  such  as  LiAl  and  NaAl  have  singlet  ground 
electronic  states  (see  below).  Finally,  our  calculated  dissociation  energy  (De)  of 
NaB  (^rir)  is  0.76  eV  at  the  QCISD(T)/6-31l+G(2df)  level. 

NaC.  For  this  molecule,  one  might  anticipate  any  of  three  valence  orbital 
occupancies  lo-  ^c^2a'^^K^  and  When  the  In  orbital  is 

occupied  by  three  electrons  (essentially  2p  on  C).  we  obtain  a  state.  When 
two  electrons  occupy  the  1n  orbital  and  the  third  electron  occupies  the 
antibonding  2o  orbital  (essentially  2s2p  hybrid  orbital  on  C).  a  state  results. 
Finally,  when  two  electrons  occupy  the  2a  orbital  and  one  electron  occupies  the 
1n  orbital,  a  ^Ur  state  results.  Among  all  states  with  n^  and  n^  occupancies, 
we  found  the  hioh-spin  (1a2ln22o^)  state  to  be  the  lowest.  The  low-spin  ^fir 
(1a22a2lni)  and  2ni  (1o2ln3)  states  are  less  stable  by  29.4  kcal/mol  and  44.8 
kcal/mol,  respectively  (see  Table  111).  We  feel  confident  in  predicting  that  is 
the  ground  electronic  state  of  NaC.  Because 

(Table  III  near  here  please) 

(i)  the  results  for  relative  energies  of  these  three  electronic  states  agree  well 
each  other  at  the  QCISD(T)  and  PMP4  levels,  (ii)  the  spin  contaminations  in  all 
three  states  are  low;  and  (iii)  the  energy  of  the  lowest  excited  2nr  electronic 
state  is  more  than  1  eV  (the  accuracy  of  relative  energies  obtained  in  our 
calculations  at  the  QCISD(T)/6-31 1+G(2df)  level  is  ca.  0.3  eV).  The  valence 
isoelectronic  LiC,  LiSi  and  NaSi  diatomics  also  have  ('*1')  ground  electronic 
states.  Our  calculated  dissociation  energy  (De)  of  NaC  is  1.97  eV  at  the 
QCISD(T)/6-31 1  +G(2df)  level. 


NaN.  For  this  molecule,  three  low-lying  states  have  been  studied; 
(1o22o2l7t2).  3^1  (icj22oiijt3)  and  ^1+  {lo2ijt4).  (Table  IV). 

(Table  IV  near  here  please) 

At  the  PMP4  and  QCISD(T)  levels,  the  is  well  separated  from  the  others  in 
energy.  Because  the  state  is  lower  in  energy  than  the  state  at  the  PMP4 
level  while  the  ^1:  state  is  lower  than  the  state  at  the  QCISD(T)  level,  we 
carried  out  large  scale  CASSCF-MRCISD(T)  calculations  using 
(17s12p5d4f/7s5p3d2f)Na  +  (14s9p4d3d/7s5p3d2f)N  basis  sets  for  these  two 
states.  The  CASSCF  expansions  included  378  (^X*)  and  384  (^Oi ) 
configurations,  which  represent  all  possible  occupations  of  six  valence 
electrons  in  (4,2,2, 1)  active  orbitals.  All  single-  and  double-excitations  from 
these  378  and  384  CASSCF  configurations  were  then  included  in  the  MRCISD 
calculations  giving  1,115,606  (^X*)  and  1,154.220  (^Ili )  MRCISD 
configurations.  Davidson's  corrections  to  the  MRCISD  energies  were  very  small 
and  thus  our  results  at  the  CASSCF-MRCISD  and  CASSCF-MRCISD(Q)  levels 
are  very  close  each  other  (Table  IV). 

The  state  was  predicted  to  be  the  ground  state  at  both  levels  with  an 
excitation  energy  of  ca.  8  kcal/mol  (see  Table  IV).  Both  the  and 

the  states  are  well  represented  by  Hartree-Fock  wave  functions 
(CHF(1<722o2l7tx“1Jty2)=0-9568  for  and  CHF(1o22o2lKx«lKya)=0.9452  for 
3X')-  From  more  than  1 ,000,000  configurations,  only  the  Hartree-Fock 
configurations  have  coefficients  larger  than  0.15.  The  occupancies  of  the 
natural  MO,  for  the  ^Fli  and  ^X'  states  are:  io2  O02a^  ^6300.99400.02 

27:x0  0l27ty0  05  and  1o2  OO20l  .963o1.8940O.O8trt^O.99i;CyO-99 
27tx9  9^2ny9-9^  respectively. 


9 


Because  the  state  has  the  lowest  energy  relative  to  the  other  states 
using  all  our  sophisticated  ab  initio  methods  (QCISD(T),  CASSCF-MRCISD  and 
CASSCF-MRCiSD(Q)V  we  are  confident  that  the  state  is  the  ground 
electronic  state  of  MaN.  The  valence  isoelectronic  LiN  molecule  also  has  a 
ground  state. Our  calculated  dissociation  energy  (De)  of  NaN  (^v-)  jg  q  77 
eV  at  the  QCISD(T)/6-31 1+G(2df)  level. 

MgB.  This  molecule  is  isoelectronic  with  BeB  which  we  discussed 
earlier.  Based  on  our  BeB  findings,  we  studied  the  following  electronic  states: 
2nr(1cT22a2l7ri),  2S+  {^a^2a^3a^).  ^flr  (1o22oi1ni3ai),  (1o22oiU2)  and 

^rii  (1o2i7t3).  Preliminary  calculations  on  these  states  at  the  MP2(full)/6*31 1+G* 
level  predict  the  ^Ur  state  to  be  the  lowest  (Table  V).  The  2nr  state  was  also 
found  to  be  the  ground  state  when  larger 

(Table  V  near  here  please) 

basis  sets  and  more  sophisticated  correlation  methods  (  PMP4  and  QCISD(T) 
with  6>31 1+G(2df)  basis  sets)  were  used.  The  2X*  state  is  predicted  to  be  the 
first  excited  state,  with  the  ^Flr  and  states  lying  higher  in  energy. 

For  MgB  we  are  confident  that  the  2nr  state  is  the  ground  electronic 
state  because  (i)  all  four  low-lying  states  have  small  spin-contamination,  (ii)  the 
energy  difference  between  the  ground  state  and  the  first  excited  state  is  7-9 
kcal/mol,  and  (iii)  the  relative  energy  difference  between  results  at  the  PMP4 
and  QCISD(T)  levels  is  only  2.4  kcal/mol,  which  is  less  than  the  first  excitation 
energy.  Our  calculated  dissociation  energy  (De)  of  MgB  (2nr)  is  0.47  eV  at  the 
QCISD(T)/6-3ll+G(2df)  level. 
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MgC.  Given  a  total  of  six  valence  electrons  (3s2  from  Mg  and  2s22p2 
from  C).  the  most  likely  candidates  for  low-energy  configurations  involve 
1o22a2ln2,  1a22oMn:2.  or  1  a22oi  1 7t23o^  occupancies.  We  therefore 

studied  the  following  five  low  lying  electronic  states  for  MgC:  (Ia22cr2ir:2) 

^rii  (10220^1 7:3),  1V+  (1021^4),  Sfir  (1o22a2l  jti3oi)  and  ^1' 
and  the  results  of  our  calculations  appear  in  Table  VI. 

(Table  Vi  near  here  please) 

We  find  the  (lo22o2irt2)  state  to  be  the  ground  state  and  the 
(1(j22oM7t23o^)  state  to  be  the  lowest  lying  excited  state  at  both  the  PMP4  and 
QC1SD{T)  levels.  However,  the  adiabatic  excitation  31-  energy  varies 
from  1.0  kcal/mol  at  the  PMP4  level  to  10.5  kcal/mol  at  the  QCISD(T)  level. 
Therefore  we  also  carried  out  calculations  on  these  two  states  at  the  CASSCF- 
MRCISD(Q)  level  using  a  (17s12p5d4f/7s5p3d2f)Mg  +  (14s9p4d3d/7s5p3d2f)c 
basis  set.  The  CASSCF  expansion  included  260  and  378  (3X*) 
configurations,  which  represent  all  possible  distribution  of  the  six  valence 
electrons,  among  the  (4,2,2, 1)  and  (3,2,2, 1)  active  orbitals,  respectively.  All 
single*  and  double-excitations  from  these  260  and  378  CASSCF  configurations 
were  then  included  in  the  MRCISD  calculations  giving  the  852,680  (^I')  and 
1,155,606  (3X-)  MRCISD  configurations.  Davidson's  corrections  to  the  MRCISD 
energy  were  modest  so  results  at  the  CASSCF-MRCISD  and  CASSCF- 
MRCISD(Q)  levels  are  reasonably  close  (Table  VI).  Because  the  3X'  state  has 
the  lowest  energy  at  all  four  of  our  most  sophisticated  ab  initio  levels,  we  are 
confident  that  3£*  is  the  ground  electronic  state  for  MgC.  Our  calculated 
dissociation  energy  (De)  for  MgC  (^L*)  'S  1.50  eV  at  the  QCISD(T)/6-31 1-«-G(2df) 


level. 


After  our  calculations  were  completed,  an  article  by  Bauschlicher, 
Langhoff  and  Partridge  (BLP)^®  appeared  in  the  literature  reporting  high  quality 
calculations  on  low-lying  electronic  states  of  MgC.  These  workers  also  found 
to  be  the  ground  state  and  to  be  the  first  excited  state.  Our  results  for  the 
ground  electronic  state  (Re(Mg-C)=2.099  A,  De=1.50  eV  and  AGi/2=527  cm- 
1)  agree  well  with  the  BLP  data  (Re(Mg*C)=2.103  A,  De=1.52  eV  and  aGi/2=541 
cm-'')  however,  for  the  first  excited  state,  our  data  ^1.'  (Re(Mg-C)=2.066  A,  De= 
eV,  Te=3975  cm*''  and  AGi/2=570  cm-1)  are  somewhat  different  from  those  of 
BLP  {Re(Mg-C)=2.109  A,  De=0.66  eV.  Te=3545  cm-''  and  AGi/2=515  cm-1).  The 
quite  substantial  difference  in  the  bond  length  (0.043  A)  and  AG  1/2  (55  cm-'') 
values  is  strange  because,  for  the  ^X-  state  we  have  good  agreement  for  these 
quantities  using  four  MP2(full),  QCISD(T),  CASSCF-MRCISD  and  CASSCF- 
MRCISD(Q)  different  methods. 

MgN.  Given  a  total  of  seven  valence  electrons  (3s2  from  Mg  and  2s^2p^ 
from  N),  the  most  likely  candidates  for  low-energy  configurations  involve 
1o22cT2i7t23<ji^  1 0^20^1  1  <y^2a^  1  and  1o^2a^lK'^  occupancies.  We 

therefore  studied  the  following  four  low-lying  electronic  states; 

(1o22o2i^23cj1).  2ni  (1a22o2ln3).  4ni  (1a22ail7i33oi)  and  2X+  (1c22(jUn^). 

The  results  of  our  calculations  appear  in  Table  VII. 

(Table  VII  near  here  please) 

We  found  the  high-soin  (1o22o2 17:23a'')  state  to  be  the  ground 
electronic  state  and  the  20)  (la22o2l7:3)  state  to  be  the  lowest  lying  excited 
state.  However,  the  adiabatic  X'^X*  ^flj  excitation  energy  is  not  large  and 
varies  from  1 1 .2  kcal/mol  at  the  PMP4  level  to  6.3  kcal/mol  at  the  QCISD(T) 
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level.  Therefore,  we  calculated  the  two  states  at  the  CASSCF-MRCISD(Q)  level 
using  a  (17s12p5d4f/7s5p3d2f)Na  +  (14s9p4d3d/7s5p3d2f)B  basis  set.  The 
CASSCF  expansion  included  104  and  196  (^Fli)  configurations,  which 
represent  all  possible  distribution  of  the  7  valence  electrons  (two  electrons 
occupy  the  lo-MO  in  all  cases)  among  the  (3.2. 2,0)  active  orbitals.  All  single- 
and  double-excitations  from  the  104  (4X-)  and  196  (Sfli)  CASSCF 
configurations  were  then  included  in  our  MRCISD  calculations  giving  955,575 
{^1.')  and  1 ,204,664  (20;)  MRCISD  configurations. 

The  results  are  presented  in  Table  VII.  The  state  has  a  lower  energy 
than  the  Sfli  state  at  both  the  CASSCF-MRCISD  and  CASSCF-MRCISD(Q) 
levels.  The  and  2ni  states  are  not  well  represented  by  the  Hartree-Fock 
wave  function  CHF(1o^2o23o''1jtx“1?Cy“)=0.3137  and 
C(1a22ai3o2l7Cx“1rty“)=0.8770  for  4$;-  and  CHF(1a22a2l7rx“l7ry2)=0.8327  for 
^rii)  in  the  MRCISD  wave  function  at  the  optimal  bond  lengths.  At  R(Mg- 
N)=2.100  A  and  shorter  distances,  the  dominant  configuration  in  the  state  is 
CHF(1<y^2cy23<Til7tx“1rcy“)=0.9526  and  all  other  configurations  in  MRCISD 
expansion  have  coefficients  less  than  0.15  in  magnitude.  At  the  equilibrium 
intern uclear  distance,  the  (1a22a23ail7tx“1ny“)  configuration  has  a  coefficient 
of  0.3137  while  the  C(1cr22ai3a2licx“1iry“)=0.8770  configuration  is  dominant 
with  all  others  having  amplitudes  less  than  0.15.  The  Davidson  correction  for 
this  state  is  also  different  at  short  and  long  interatomic  distances  as  a  result  of 
which  the  findings  at  the  CASSCF-MRCISD  and  CASSCF-MRCISD(Q)  levels 
are  quite  different. 

For  the  first  excited  2ni  state,  the  leading  configuration 
(CHF(1<y^2o23a0l7tx“1iry2)=0.8327)  remains.the  same  along  the  potential 
energy  curve,  but  two  other  configurations  C(1 0220^20^1  nx“l7ry^)=-0. 3482  and 
C(1a22o«2a“l7txPl7ty2)=-0.2806  have  amplitudes  larger  than  0.15.  The 
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occupancies  of  the  natural  MOs  for  the  and  2ni  states  are: 

1(ji  .962ai  9i3o0.994a0  03  irtxO  99i;ry0.99  and 

1  a''  962o1  603(50.384(j0.02i 7tj^o.98i 7tyi.9i  2jtx°  9^,  respectively.  Because 
the  '♦S*  state  has  the  lowest  energy  relative  to  the  other  states  at  all  four  of  our 
most  sophisticated  ab  initio  levels,  we  are  confident  that  the  state  is  the 
ground  electronic  state.  Our  calculated  dissociation  energy  (De)  of  MgN  is 
0.50  eV  at  the  QCISD(T)/6-31 1+G(2df)  level. 

AIB.  The  two  valence  isoelectronic  molecules  B2  and  AI2  are  known  to 
have  and  ground  electronic  states,  respectively.  However,  the 

B2  molecule  has  a  very  low-lying  ^Xu'  first  excited  state  (Te=1701  cm'TSi)  and 
AI2  has  a  low-lying  ^Xg*  first  excited  state  (Te=200  cm*T32)_  Therefore,  for  AIB 
we  anticipate  that  one  of  the  ^Or  (1cr22o2l7r''3o‘'),  ^X*  (1o22o'' 1^230“')  or  3X' 
(1cj22o21k2)  states  will  be  the  ground  state,  while  the  others  are  low  lying 
excited  states.  We  carried  out  calculations  on  states  of  these  three  symmetries 
at  several  levels  of  theory  (Table  VIII). 

(Table  VIII  near  here  please) 

At  the  MP2(full)/6-31 1+G*  level,  the  ^Or  state  is  the  lowest  with  the  ^X' 
the  first  and  the  ^X'  the  second  excited  states.  We  find  the  same  ordering  at  the 
PMP4  and  QCISD(T)  levels;  however,  at  the  QCISD(T)/6-31 1+G(2df)  level,  the 
energy  difference  between  the  ^Or  and  ^X'  states  is  only  0.2  kcal/mol.  The  ®X' 
second  excited  state  is  higher  in  energy  by  22.6  kcal/mol,  and  therefore,  is  not  a 
candidate  for  the  ground  electronic  state  of  AIB.  Although  the  state  is  lowest 
at  all  levels  of  theory,  the  spin-contamination  of  the  ^X'  state  is  very  high. 
Because  these  states  have  very  nearly  identical  total  energies  at  both  the  PMP4 


and  QCISD(T)  levels,  we  are  not  able  to  predict  with  certainty  the  ground 
electronic  state  from  these  data.  Therefore,  we  studied  these  states  again  at  the 
MCSCF  and  CASSCF-MRCISD(Q)  levels  of  theory. 

MCSCF  calculations  have  been  performed  using  6  valence  electrons  in 
8  valence  active  MOs  (giving  1512  CSFs).  The  three  lowest  triplet  roots  were 
then  calculated  using  the  SA  MCSCF  methodology.  The  optimized  bond 
lengths  for  the  ^Or  and  states  at  this  level  are  very  close  to  the  those  at  the 
QCISD(T)  level,  while  the  MCSCF  harmonic  frequencies  are  both  lower  than  at 
the  MP2(full)/6-31 1+G*  level.  The  state  is  the  lowest  at  MCSCF  level,  and 
the  is  the  first  excited  lying  2.6  kcal/mol  higher.  This  ordering  contradicts  all 
of  our  previous  results.  Although  these  MCSCF  calculations  have  no  spin- 
contamination  and  all  important  valence  orbitals  are  included  in  variational 
calculation,  the  fraction  of  the  total  correlation  energy  included  at  the  MCSCF 
level  is  not  high.  Therefore,  we  examined  these  two  states  at  the  CASSCF- 
MRCISD(Q)  level  using  a  (17s12p5d4f/7s5p3d2f)Ai+  (14s9p4d3d/7s5p3d2f)B 
basis  set.  The  CASSCF  expansion  included  378  and  384  (SHr) 
configurations,  which  represent  all  possible  distribution  of  6  valence  electrons 
in  (3,2,2, 1)  active  orbitals.  All  single-  and  double-excitations  from  these  378  (3t- 
)  and  384  (^Flr)  CASSCF  configurations  were  then  included  in  the  MRCISD 
calculations  giving  1,155,606  (^S*)  and  1,154,220  (^Or)  MRCISD 
configurations.  Results  of  these  calculations  are  presented  in  Table  VIII. 

The  31'  state  is  predicted  to  be  most  stable  at  both  the  CASSCF- 
MRCISD  and  CASSCF-MRCISD(Q)  levels.  The  ->  ^Or  excitation  energy  is 
2.4  kcal/mol  (CASSCF-MRCISD)  and  1.6  kcal/mol  (CASSCF-MRCISD(Q)). 
which  agree  well  with  the  MCSCF  results.  Both  the  and  ^Or  states  are  well 
represented  by  Hartree-Fock  wave  functions  CHF(1<y^2a23oOl7tx“f5^y“)=0-8952 
for  3S'  and  CHF(f  <y^2a23<j“lJCx“)=0.9061)  for  Sflr  at  the  optimal  bond  lengths. 
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From  more  than  1 ,000,000  configurations  in  MRCISD  expansions,  only  one 
non-Hartree-Fock  configuration  has  a  coefficient  (C{1o22a03o2ljrx^  ■'rry^)=- 
0.1985  for  3X‘)  larger  than  0.15.  Moreover,  the  Davidson  correction  is  small  and 
therefore  the  results  at  the  CASSCF-MRCISD  and  CASSCF-MRCISD(Q)  levels 
are  similar.  The  occupancies  of  the  natural  MOs  for  the  and  ^Flr  states  are: 
1(ji  .922ai  7®3o0  164(j0,02  ®®1rcy0-96  2rtx®  and 

1  ai  922ai  853(j0.994(j0.03i  ;jj^o.97-i  062jtx0  062ny0  03^  respectively. 

We  are  confident  that  ^1.'  is  the  ground  electronic  state  for  AIB  be 
(i)  the  results  are  the  same  at  the  MCSCF  and  CASSCF-MRCISD(Q)levels 
theory  and  (ii)  B2  has  a  ^Xg-  ground  electronic  state  with  a  1701  crn'T^i 
excitation  energy  into  ^Hu  state,  (iii)  in  AI2,  the  state  is  located  only  200 
cm-1  above  the  state. ^2  Therefore  we  expect  that  in  AIB  the  energy 
difference  between  these  two  states  should  be  somewhere  inbetween.  Our 
calculated  dissociation  energy  (Dg)  of  AIB  (^I*)  is  1.78  eV  at  the  QCISD(T)/6- 
311+G(2df)  level. 

NaAI.  In  previous  work  on  LiB,  NaB  and  LiAl,  three  states  of  symmetries 
■>1+  (1<T22a2),  3X'  {1a2lrr2)  and  ^flr  (1a2l;ri2oi)  were  found  to  be  low-lying.  In 
LiB  and  NaB,  is  the  ground  state  while,  for  the  valence  isoelectronic  LiAl, 
the  state  is  the  lowest.  We  examined  these  three  low-lying  states  for  NaAI. 
As  in  the  LiAl  molecule,  we  found  that  the  low-spin  i^-*-  state  is  the  most  stable 
for  NaAI.  However,  the  high-spin  ^fir  state  is  only  5.2  kcal/mol  above  the  ground 
state  at  the  QCISD(T)/6-31 1+G(2df)  level,  and  the  state  lies  above  the  ^1+ 
state  by  22.4  kcal/mol.  The  dissociation  energy  (Dg)  of  NaAI  CX'*’)  is  predicted 
to  be  17.8  kcal/mol  at  the  QCISD(T)/6-311 +G(2df)  level.  The  energy  difference 
between  the  ‘'X+  and  ^Flr  states  for  valence  isoelectronic  LiAl  was  studied  at  the 
CASSCF-MRCISD(Q)  level  using  a  {17s12p5d4f/7s5p3d2f)Ai  + 
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(14s9p4d3cl/7s5p3d2f)Lj  basis  set  (see  below),  and  we  obtained  a  very  good 
agreement  between  QCISD(T)  and  CASSCF-MRCISD{Q)  level  results  for  the 
^Or  excitation  energy.  Because  the  energy  differences  between  the 

and  states  for  LiAl  and  NaAl  are  the  same,  we  are  confident  that  our 
prediction  of  the  ground  electronic  state  for  NaAl  is  reliable. 

MgAI.  This  molecule  is  valence  isoelectronic  to  BeB,  so  we  studied  the 
same  five  low-lying  electronic  states:  2X+  (1a22a23oi), 

(1cT22o''1rt''3o‘'),  (1o22afl7:2)  and  2ni  {1o2l;:3)  as  we  identified  for  BeB. 

Preliminary  calculations  at  the  MP2(full)/6-311+G*  level  were  carried  out  for 
^Or.  “^rir.  ^S'.  and  ^Rj  states.  We  found  the  state  to  be  the  lowest 
(Table  X). 

(Table  X  near  here  please) 

Using  more  sophisticated  correlation  method  (QCISD(T)/6-311+G(2df)),  the  2X+ 
first  excited  state  is  found  to  be  only  4.5  kcal/mol  above  the  ground  ^R^  state. 
Therefore  we  studied  these  two  lowest  states  at  the  MCSCF  level  of  theory. 

The  MCSCF  calculations  were  performed  using  6  valence  electrons  in  8 
valence  MOs  which  gave  rise  to  1512  CSFs.  The  three  lowest  doublet  roots 
were  calculated  using  the  SA  MCSCF  methodology.  The  ^Ur  state  was  also 
found  to  be  lowest  at  the  MCSCF  level.  The  2X+  state  is  higher  in  energy  by 
15.0  kcal/mol.  Thus,  for  MgAI  we  are  confident  that  the  ^R^  state  is  the  ground 
electronic  state  because  (i)  the  two  low-lying  states  have  small  spin- 
contaminations,  (ii)  the  ground  electronic  staUe  is  the  same  at  the  PMP4, 
QCISD(T)  and  MCSCF  levels,  and  (iii)  the  relative  energy  range  between  PMP4 
and  QCISD(T)  is  only  0.8  kcal/mol,  less  than  the  first  excitation  energy.  Our 


calculated  dissociation  energy  (De)  of  MgAI  (Sflr)  is  0.34  eV  at  the  QClSD(T)/6- 
311+G(2df)  level. 


In  summary  we  have  identified  the  ground  electronic  states  of  ten 
molecules  that,  as  yet,  are  theoretically  and  experimentally  uncharacterized.  In 
the  following  section,  we  present  results  of  our  current  calculations  on  five  other 
molecules  which  were  studied  before:  LiAl,  BeP.  BP  and  SiP  but  for  which  the 
ground  electronic  states  have  not  been  determined  with  certainty. 

B.  Five  Oiatomics  Whose  Ground  States  Are  Uncertain 

LiAi.  In  our  previous  work^ioc  vve  studied  three  (1o22a2),  3£' 
(1o2ln2)  and  ^rir  (1o2l7t‘'2a'')  low-lying  states  of  LiAl  (Table  XI),  and 

(Table  XI  near  here  please) 

found  the  low-spin  state  to  be  the  most  stable.  However,  the  lowest  high- 
spin  ^rir  state  is  predicted  to  be  only  5.1  kcal/mol  higher  at  the  QCISD(T)/6- 
31 1+G*  level,  and  the  state  lies  only  21.2  kcal/mol  above  the  Our 
calculated  dissociation  energy  for  LiAl  (iZ*^)  was  23.3  kcal/mol. 

Because  the  energy  difference  between  the  and  ^flr  states  is  small, 
we  performed  calculations  on  these  two  states  at  the  CASSCF-MRCISD(Q) 
level  using  (17s12p5d4f/7s5p3d2f)Na  +  (14s9p4d3d/7s5p3d2f)B  basis  sets.  The 
CASSCF  expansion  included  152  (^Z"^)  and  160  (^Flr)  configurations,  which 
represent  all  possible  distribution  of  the  valence  electrons  among  (4,2,2, 1) 
active  orbitals.  All  single-  and  double-excitations  from  these  CASSCF 
configurations  were  then  included  in  our  MRCISD  calculations  giving  84,149 
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snd  117,475  (^Hr)  MRCISD  configurations.  The  Davidson  correction  to 
the  MRCISD  energy  was  very  small  for  both  states,  so  our  results  at  the 
CASSCF-MRCISD  and  CASSCF-MRCISD(Q)  levels  are  very  similar  (Table  XI). 

The  state  is  found  to  be  the  most  stable  at  both  the  CASSCF- 
MRCISD  and  CASSCF-MRCISD(Q)  levels,  with  the  ^Flr  state  higher  in  energy 
by  5.0  kcal/mol.  This  is  the  same  excitation  energy  we  found  at  the  QCISD(T) 
level,  so  we  are  confident  that  ‘•1+  is  the  ground  electronic  state  for  UAL  Both 
the  ^Or  and  “'2,+  states  are  well  represented  by  Hartree-Fock  wave  functions 
(CHF(1<y^2o2)=0.9163  for  12+  and  CHF(1o^2a'' l7tx“)=0.9476  for  Sflr)-  Among  all 
configurations  in  the  MRCISD  expansions  for  both  these  configurations,  only 
one  configuration  has  a  coefficient  (C(1a22o03o2)=-0.1656  for  32*)  larger  than 
0.15.  The  occupancies  of  the  natural  MO  for  the  ■'2‘^  and  ^Y\r  states  are; 

1(yi.882<ji  783(j0.08i jjj^o.osi TCyO.oa  27tx®'®^2)iy0 03  and 

1(ji.882a0.983<j0.02i7Cj^0.97i;fy0.05  2TtxO  032;ty0.oi  ^  respectively. 

MgSi.  In  our  previous  work^it)  vve  studied  five  ■'2'*’  (la^lrr^),  3£- 
{1a22o2l7i2),  3nr  i^a22c2^K^3<s^),  3ni  (1o22ail7i3)  and  {la22ail7t23oi) 
electronic  states  of  the  MgSi  molecule,  and  we  found  the  3^*  state  to  be  the 
most  stable.  However,  the  high-spin  3nr  state  is  only  2.7  kcal/mol  less  stable  at 
the  PMP4/6-31 1+G*  level.  Because  the  energy  difference  between  the  3X'  and 
3nr  states  is  small,  we  carried  out  calculations  on  these  two  states  at  the 
CASSCF-MRCISD(Q)  level  using  {17s12p5d4f/7s5p3d2f)Mg,Si  basis  sets.  The 
CASSCF  expansion  included  378  (32*)  and  384  (3nr)  configurations,  which 
represent  all  possible  distribution  of  the  valence  electron  among  (3,2,2, 1)  active 
orbitals.  All  single-  and  double-excitations  from  these  CASSCF  configurations 
were  than  included  in  the  MRCISD  calculations  giving  1,054,150  (32*)  and 
1 ,052,968  (^Or)  MRCISD  configurations.  The  Davidson  correction  to  the 


MRCISD  energy  was  very  small  for  both  states  as  a  result  of  which,  results  at 
the  CASSCF-MRCISD  and  CASSCF-MRCISD(Q)  levels  are  very  close  to  each 
other  (Table  XII). 

The  state  is  predicted  to  be  the  most  stable  at 

(Table  Xll  near  here  please) 

both  the  CASSCF-MRCISD  and  CASSCF-MRCISD{Q)  levels,  and  the  state 
is  higher  in  energy  by  16-18  kcal/mol.  Therefore  we  are  confident  in  suggesting 
that  is  the  ground  electronic  state  for  MgSI.  Both  the  ^Flr  and  states  are 
quite  well  represented  by  a  Hartree-Fock  wave  function 
(CHF(1o22o23af  l7tx^)=0.9035  for  and  CHF(1o22a2ljtx«l7Cy«)=0.8851  for 
aZ').  Among  all  of  the  configui  itions  in  the  MRCISD  expansions  for  both  these 
configurations,  the  only  configuration  with  a  coefficient  larger  than  0.15  is 
(C(1o22oQ3o2ijrx“lTty“)=-0.2459  for  a^-).  The  occupancies  of  the  natural  MO 
for  the  and  aplr  states  are:  lo^  a22aiT53(j0.2i4^y0.02i;ij^0.97i;jy0.97 

2KxO  oa2jty^  oa  and  1  oi  31201  86300.98400.04-1  ;jj^0.98i7jy0.07  27:x®  ®^27tyO  oa 

respectively. 

BeP.  In  our  previous  study^4  w©  found  two  electronic  states 
(1<T22a2l7c23a^)  and  2ni  (1a22a2lna}  to  be  the  most  stable  for  BeP.  At  the 
PMP4/6-31 1-HG(2df)  and  at  QCISD(T)/6-31 1+G(2df)  levels,  the  ^1*  state  is 
lower  by  only  3. land  3.5  kcal/mol,  respectively.  Calculations  on  these  two 
states  at  the  CASSCF-MRCISD(Q)  level  using  a  (17s12p5d4f/7s5p3d2f)p + 
{14s9p4d3d/7s5p3d2f)Be  basis  set  produced  expansions  with  104  (4X‘)  and 
196  (2ni)  configurations,  which  represent  all  possible  distributions  of  7  valence 
electrons  in  (3,2,2,0)  active  orbitals.  All  single-  and  double-excitations  from 
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these  CASSCF  configurations  were  then  included  in  MRCISD  calculations 
giving  955,575  and  1,204,664  (Sflj)  configurations.  The  Davidson 
correction  to  the  MRCISD  energy  was  modest  for  both  states,  so  results  at  the 
CASSCF-MRC.3D  and  CASSCF-MRCISD(Q)  levels  are  reasonably  similar 
(Table  Xlll). 

The  4X*  state  is  found  to  be 

(Table  Xlll  near  here  please) 

the  most  stable  at  both  the  CASSCF-MRCISD  and  CASSCF'MRCISD(Q) 
levels.  Because  this  state  is  the  ground  electronic  state  at  our  four  most 
sophisticated  ab  initio  methods  (PMP4,  QCISD(T),  CASSCF-MRCISD  and 
CASSCF-MRCISD(Q)),  we  are  confident  that  4X*  is  the  ground  electronic  state 
for  BeP. 

Both  the  '*I‘  and  2nj  states  are  well  represented  by  Hartree-Fock  wave 
functions  (CHF(1cr^2a23cr«l7tx“1rty®)=0.9437  for  4X*  and 
CHF(1<y^2a2l7Cx“1rty2)=0.9098  for  200.  Among  all  of  the  configurations  in  the 
MRCISD  expansions  for  both  these  states,  only  the  Hartree-Fock  configurations 
have  coefficients  larger  than  0.15.  The  occupancies  of  the  natural  MO  for  the 
4^*  and  20]  states  are:  la''  ^23^0.994(70.031  71^0,99  2nx®-®^2jty0 

and  1a^  ®52aT853(j0.i04(j0.02i7tx0  97ijryi.88  27rx° °'*27tyO 08^  respectively. 


BP.  Previously,'' 4  we  found  two  low  energy  electronic  states 
(1a22a2l7t4)  and  ^l\\  for  BP^  The  ^Fli  state  is  lower  by  only  1.8 

and  6.8  kcal/mol,  at  the  PMP4/6-31 1+G(2df)  and  at  QCISD(T)/6-31 1+G(2df) 
levels,  respectively.  We  therefore  carried  out  calculations  of  these  two  states  at 
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the  CASSCF-MRCISD(Q)  level  using  (17sl2p5d4f/7s5p3d2f)p  + 
(14s9p4d3d/7s5p3d2f)Be  basis  sets.  The  CASSCF  expansion  included  142 
CS'*’)  and  150  (^Oi)  configurations,  which  represent  all  possible  distribution  of 
the  7  valence  electrons  (two  electrons  occupy  the  1a-MO  in  all  cases)  among 
the  (3, 2, 2,0)  active  orbitals.  All  single-  and  double-excitations  from  these 
CASSCF  configurations  were  then  included  in  our  MRCISD  calculations  giving 
1,089,382  CS"^)  and  1,775,350  (^Fli)  configurations.  The  Davidson  correction  to 
the  MRCISD  energy  was  modest  as  a  result  of  which  the  CASSCF-MRCISD 
and  CASSCF-MRCISD(Q)  results  are  similar  (Table  XIV)  but  not  identical. 

(Table  XIV  near  here  please) 

The  ^rii  state  is  predicted  to  be  the  most  stable  at  both  the  CASSCF- 
MRCISD  and  CASSCF-MRCISD(Q)  levels.  Because  the  ^Fli  state  is  the  lowest 
at  all  sophisticated  ab  initio  levels  (PMP4,  QCISD(T),  CASSCF-MRCISD  and 
CASSCF-MRCISD(Q))  we  are  reasonably  confident  that  is  the  ground  state 
for  BP. 

Both  the  and  ^Oi  states  are  well  represented  by  Hartree-Fock  wave 
functions  (CHF(1<J^2o2i7c^2i7jy2)=o.8768  for  1^+  and 
CHF{1o^2o22oil7:x“1«y^)=  0.9123  for  Among  all  configurations  in  the 
MRCISD  expansions  for  both  of  these  states,  only  one  non-Hartree-Fock 
configuration  (C(1a22o03o2lKx2l7ty2)=-0.1654  for  1I-")  has  a  coefficient  larger 
than  0.15.  The  occupancies  of  the  natural  MOs  for  the  and  states  are; 

1  oi  .962o1  823cy0. 1  64<j0.02  1  t^^i  .86i  86  ^  ^  27tyO  ^  ^  and 

1  oi  962o1  .903<j1  .004(y0.03i  kx^-97 i  ^^1  8727t^o.o627ty0.i 0,  respectively. 
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SIP.  In  our  previous  study"''*  we  identified  two  low  energy  electronic 
states  (1o22a2l;t'*3o')  and  2ni  {^a^2a^^n^Zo^)  for  the  SiP  molecule.  At 
the  PMP4/6-31l+G(2df)  level,  the  state  is  lower  by  3.3  kcal/mol,  but  at  the 
QCISD(T)/6-31 1-*-G(2df)  level,  the  2nj  is  more  stable  by  0.9  kcal/mol.  Again,  we 
carried  out  calculations  of  these  two  states  at  the  CASSCF-MRCISD(Q)  level 
using  (17s12p5d4f/7s5p3d2f)si,P  basis  sets.  The  CASSCF  expansion  included 
208  (25;+)  and  196  (2nj)  configurations,  which  represent  all  possible  distribution 
of  the  7  valence  electrons  (two  electrons  occupy  the  I0-MO  in  all  cases)  among 
the  (3,2.2, 0)  active  orbitals.  All  single-  and  double-excitations  from  these 
CASSCF  configurations  were  then  included  in  our  MRCISD  calculations  giving 
420,340  (2X+)  and  416,604  (2ni)  configurations.  The  Davidson  correction  to  the 
MRCISD  energy  was  modest  so  the  CASSCF-MRCISD  and  CASSCF- 
MRCISD(Q)  results  are  in  reasonably  close  agreement  (Table  XV). 

(Table  XV  near  here  please) 

The  2ni  state  is  found  to  be  the  most  stable  at  both  the  CASSCF- 
MRCISD  and  CASSCF-MRCISD(Q)  levels.  Because  the  201  ground  electronic 
state  is  lowest  at  three  sophisticated  ab  initio  methods  (QCISD(T),  CASSCF- 
MRCISD  and  CASSCF-MRCISD(Q))  we  are  reasonably  confident  in  claiming 
that  is  the  ground  electronic  state  for  SIP.  Both  the  2^+  and  2ni  states  are 
well  represented  by  Hartree-Fock  wave  functions 
(CHF(1o22o23cTi1jtx2ljry2)=0.9042  for  ^1*  and 

CHF(1<J^2o23o2l7tx“l7ry2)=0.9187  for  20, ).  Among  all  configurations  in  the 
MRCISD  expansions  for  both  of  these  statesr  only  the  Hartree-Fock 
configurations  have  coefficients  larger  than  0.15.  The  occupancies  of  the 
natural  MOs  for  the  ^1*  and  2ni  states  are: 
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1o2  002a^  903o1  004a0.03i;j^1.87i;jy1.87  23txO  102jty0  ^0  and 
1  o2.002o1  943q1  .904<j0.97i  ,ijj0.97i  jtyi.87  2!tx0  0®2jty0  ^0,  respectively. 

IV.  OVERVIEW 

In  Figure  2.  the  term  symbols  of  the  ground  electronic  states  calculated 
and  experimentally  determined  (where  known),  are  given  for  all  120  first  and 
second  row  diatomic  molecules,  and  in  Table  XVI  the  theoretical  and  available 
experimental  dissociation  energies  for  these  ground  states  are  presented. 

(Figure  2  and  Table  XV I  near  here  please) 

In  Table  XVi  some  of  the  experimental  dissociation  energies  are  De  values  and 
some  are  Do.  Because  the  differences  between  such  values  are  smaller  in  most 
cases  than  the  accuracy  of  our  calculations,  we  do  not  emphasize  these 
differences.  Moreover,  we  do  not  site  ail  of  the  theoretical  data  published  in  the 
literature  but  select  what  we  feel  is  the  most  reliable  data. 

The  conventional  valence  MO  ordering  for  first-row  homonuclear 
diatomic  molecules  is:  10g<1<Tu<l7tu<2og<lJCg<2au  lithium  to  nitrogen,  with 
the  Ijiu  and  2ag  orbitals  reversed  for  O2  and  Fz.^  From  this  ordering  one  can 
successfully  predict  the  ground  electronic  states  for  all  first  and  second  row 
homonuclear  diatomics,  and  these  predictions  agree  with  the  findings  given  in 
Fig.  2. 

However,  when  heteronuclear  species  are  considered,  it  is  difficult  to 
predict  the  order  in  which  the  valence  MOs  are  filled  especially  for  the  4-8 
valence  electron  cases.  Assuming  an  MO  ordering  analogous  to  that  found  in 
most  homonuclear  cases:  1o<2o<l7c<3o<2it<4o,  one  would  predict  the 
following  ground  electronic  states  for  2-16  valence  electrons  :  (1<7^). 
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(1o22o<).  (1o22o2).  mr  (1o22o2ljti).  (lo22a2l jc2),  ajli  (1<522o2u3), 

11+  (1o22a2u^),  25:+  (la22o2l7c^3oi).  1I+  (1a22a2lji43o2),  20^ 
(1a22a2l7t43<j227ti),  35;-  (Ia22a2ljc43a22;c2),  2r]i  (1a22a2l7t43cf22jc3)_  15;+ 

(1o22a2l;i43<j22ji4)^  25;^+  (Ia22a2l  Ji43a227t44ai)  and  1I+ 

(1  a22a2l  7t43a22rt44o2). 

As  shown  in  Fig.  2  in  green,  most  first  and  second  diatomic  molecules 
have  the  ground  electronic  states  that  are  expected  based  on  the  above  MO 
ordering; 

H2.  LiH.  Lia.  LiNa,  Naa  (II+  (1o2)). 

BeH,  LiBe,  MgH,  LiMg,  NaBe,  NaMg  (2l+  (1o22ai)), 

BH.  Bea,  AlH,  LiAl,  BeMg,  NaAl,  Mga  (^1+  (1c522a2)), 

CH,  BeB,  SiH,  BeAl.  MgB,  MgAI  (2nr  {1a22o2l;ci)), 

NH,  LiN,  BeC,  Ba.  PH.  UP.  BeSi,  BAl.  MgC,  NaN,  NaP,  MgSi  {^I* 
(1(j22o2u2)),  oh,  LiO,  SH,  US.  NaO,  NaS  (201  (1o22o2ln3)), 

HF.  LiF,  BeO,  Ca,  HCI.  LiCI,  BeS,  MgO.  NaF,  NaCI,  MgS  (II+ 

(1o22o21«4)),  BeF,  BO.  CN.  BeCI,  BS.  PC,  SiN,  AlO,  MgF  (21+ 

(1a22o2iji43<yi))^ 

BF,  CO,  Na,  BCI,  CS.  PN,  SiO,  AIF  CIAI,  SiS,  Pa  (^1+  (1o22a2l7i43a2)), 
CF,  NO.  CCI,  SN.  PO.  SiF,  SiCI.  SP  {2nr  (1o22a2l 71^3022711)), 

NF.  O2.  NCI,  SO.  PF,  PCI,  Sa  (3X-  (1a22o2l7i43o227c2)), 

OF.  CIO,  SF,  SCI  (2ni  (1o22o2l7i43a227t3)) 
and  Fa,  CIF,  CI2  (^I-" 

However,  twenty  diatomIcs  (marked  in  red  in  Fig.  2)  do  not  fit  the  pattern  and 
have  unexpected  high-spin  ground  electronic  states. 

The  possibility  for  ground  states  with  higher  spin  multiplicity  than 
expected  considering  the  Aufbau  principle  arises  due  to  the  near  degeneracy  of 
the  2a-  and  1  tc-MOs  or  of  the  1  tc  -  and  3o  MOs.  For  example,  the  1  a22o2i  TtO 
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states  lie  slightly  above  the  corresponding  1o22o11ji<  states  for  LiB  and  NaB, 
and  the  Icy22a2ljtf  states  are  above  the  la22aM«2  states  for  LiC,  NaC,  LiSi, 
and  NaSi. 

Another  group  of  high-spin  ground  state  molecules  anse  because  of 
quasi-degeneracy  between  the  1jc  and  3o-orbitals.  For  example,  1o22a2l;i3 
lies  above  1o22o2ix23o^  for  BeN,  BC.  BeP,  BSi,  AlC,  MgN,  MgP,  AlSi,  and 
1o22o2l7t4  states  lie  above  1o22g2iji33<yi  for  BN,  BP,  CSi,  AIN,  AlP,  and  Si2. 
These  two  groups  of  "peculiar"  species  arise  when  there  are  4-5  or  7-8  valence 
electrons. 

Clearly,  the  isoelectronic  phnciple34,  which  is  widely  used  in  chemistry, 
does  not  work  well  for  certain  diatomic  molecules  with  4-6  valence  electrons.  It 
appears  that  the  Aufbau  principle,  using  the  "standard"  order  of  MOs  works  for 
all  first  and  second  row  diatomics  except  some  of  those  with  4-5  or  7-8  valence 
electrons.  However,  not  all  of  the  4-5  or  7-8  electron  cases  produce  high-spin 
(i.e.  2<jiljt^  2<jh7t2,  1x23o^  or  lic33o^)  ground  states;  some  of  the  4-5  electron 
cases  yield  2a2  or  2a^^n^  low-spin  ground  states,  and  some  of  the  7-8  electron 
molecules  produce  or  1]c4  ground  states.  Inspection  of  Figures  3,  where 
the  low-  and  high  spin  ground  states  are  examined  according  to  their  number  of 
valence  electrons,  reveals  that  when  group  3, 4,  and  5  elements  are  involved, 
high-spin  states  can  be  expected;  group  6  and  7  elements  do  osi  produce  high- 
spin  ground  states.  We  therefore  note  in  closing,  that  challenges  remain  for 
students  of  chemistry  even  within  the  realm  of  predicting  the  energy  ordering  for 
electronic  states  of  diatomic  molecules  and  for  understanding  why  such 
ordering  arise.  Even  greater  challenges  arise  when  consideriny  three-  or  larger 
atomic  clusters  especially  when  group  3. 4  o'r  5  elements  are  involved. 
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Figure  Captions 

Figure  1 .  The  experimentally  determined  ground  electronic  states  of 
diatomic  molecules  composed  of  first  and  second  row  atoms 
(including  hydrides).  The  diatomic  molecules  with  experimentally 
unknown  ground  electronic  states  marked  in  burgundy. 

Figure  2.  The  ground  electronic  states  of  first  and  second  row  atoms 

diatomic  molecules  based  on  experimental  and  theoretical  data. 
Diatomics  with  high-spin  ground  electronic  states  are  marked  in 
red. 

Figure  3.  Ground  states  of  first  and  second  row  diatomic  molecules  viewed 
in  terms  of  the  number  of  valence  electrons.  High-spin  exceptions 
are  noted  in  red. 


Table  I.  Calculated  Molecular  Properties  o1  the  Lowest  BeB  states. 


Table  II.  Calculated  Molecular  Properties  of  the  lowest  NaB  states. 


Table  ill.  Calculated  Molecular  Properties  of  the  lowest  NaC  states. 


NaC  (^r) 

NaC  (2n,) 

NaC  (2ni) 

1o2l  11^201 

1o22o2l7i1 

lo^lit^ 

MP2(full)/6-311+G* 


QCISD(T)/6-31UG(2df 


//QCISDfn/6-31UG(2 


=-199.74712 


MP2(fullV6-311+G* 


EmP2*- 186.56706 


0)0=443  cm’1 

0)0=401  cm*1 

0)0=474  cm'1 

<S2>=3.756 

<S2-»=0.815 

<S2>=0  799 

QCISD(T)/6-311+G(2 


//MP2(fullV6-311+G* 


QCISD('n/6-311+G(2df 


EPMP4=- 199.69560 


EQCISD('n=-1 99.69677 


<s2>=2.756 


TePMP4=8  8  kcai/mol 


=0.0  kcal/mol 


EMP4*-1 99-64374 


tQCISDfn  =**199.64996 


<S2>=0.817 


TeMP4=32.5  kcal/mol 


29.4  kcal/mol 


EmP4*-1  99.61 549 


<S2>s0.792 


TeMP4=’50.3  kcal/mol 


=44.8  kcal/mol 


Table  IV.  Calculated  Molecular  Properties  of  the  lowest  NaN  states. 


Table  V.  Calculated  molecular  properties  of  the  lowest  MgB  states. 
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Table  VI.  Calculated  Molecular  Properties  of  the  lowest  MgC  states. 


Table  VII.  Calculated  molecular  properties  of  the  lowest  MgN  states. 


Table  VIII.  Calculated  molecular  properties  of  the  lowest  AIB  states. 


A©  (3X-) 

AIB  (3nr) 

AIB  (51*) 

10^20^1 

Io22o2ijil3of 

1o^2o’ 1n23Ql 

MP2(fullV6-311+G* 

MP2(fu«)/6-311+G* 

MP2(fullV6-311-t.G* 

EmP2*-266.67110 

Em  P2=‘266. 67760 

EmP2=*266.65496 

Rs(AI-B)=2.019  A 

Ra{AI-B)=2.212  A 

R8(AI-B)=1  .943  A 

a)e=651  cm'^ 

(de=521  cm"’ 

©8=792  cm*’ 

<S^>=2.644 

<S2>a2.014 

<S2>s6.002 

MCSCF//MCSCF 

MCSCF//MCSCF 

Re(AI-B)s2.088  A 

Ra(AI-B)=2.261  A 

aje=574  cm"'' 

0)a=471  cm*’ 

EmCSCF=-266.50779 

EmcsCF=-266.50371 

TeMCSCF=0.0  kcal/mol 

TeMCSCF=2.6  kcal/mol 

QCISDI‘n/6-31 1  +G(2dfV/ 

QCISDm/6-31 1  +G(2dh// 

QCISDCn/6-31 1  +G(2dfV/ 

CX:iSD(T)/6-311+G(2df) 

QCISDfn/6-311+G(2dft 

MP2(fullV6-311+G* 

Re(AI-B)=2.056  A 

Re(AI-B)=2.226  A 

R(AI-B)=1,943  A 

EpMP4=-266.57702 

EpmP4=*266.58174 

EpMP4=*266.54922 

EQCISDfT^=-266.58794 

EqcisD(TI=*266.58820 

EqCISD(TI=*266.55221 

TePMP4^0  kcal/mol 

TaPMP4=0.0  kcal/mol 

TePMP4=’7.4  kcal/mol 

TeQClSOm=0.2  kcal/mol 

TeQClSOm=00  kcal/mol 

TeQClSDm=22.6  kcal/mol 

CASSCF-MRCISD 

CASSCF-MRCISD 

Re(AI-B)=2.046  A 

Re(AI-B)=2.224  A 

a)e=608  cm*^ 

©8=493  cm*^ 

EcaSSCF-MRCISD='266,601 94 

EcASSCF-MRCISD**266.59809 

TeCASSCF-MRCISD=0.0  kcal/mol 

TaCASSCF-MRCISD»2.4  kcal/mol 

CASSCF-MRCISD(Q) 

CASSCF-MRCISD(Q) 

Re(AI-B)=2.051  A 

Re(AI-B)=2.237  A 

(de^^OO  cm*’ 

©8=482  cm*’ 

EcASSCF-MRCISD(QI*-266.60845 

EcASSCF-MRC1SD(QI=-266.60585 

TaCASSCF.MRCISD/QI=0.0  kcal/mol 

TeCASSCF-MRCISD(QI==1  6  kcal/mol 

OV 


Table  IX.  Calculated  Molecular  Properties  of  the  lowest  NaAl  states. 


NaAl  (1l+) 

NaAl  (3nr) 

NaAl  (3r) 

1a^2a^ 

1o2l}il2o^ 

10^1^2 

MP2(fullV6-311+G* 

MP2(fullV6-3n+G* 

MP2(fullV6-311+G* 

EmP2="^04. 02294 

EmP2=-404.01938 

EmP2=-403. 98789 

Re(Na-AI)=3.134A 

Re(Na-AI)=2.985  A 

Re(Na-AI)=2.673  A 

(oe=196  cm'^ 

(0e=184  cm'^ 

c)e=259cm*^ 

<S2>=2.024 

<S2>=2.044 

QCISDCn/6-311+G(2df) 

QCISD(T)/6-311+GI2df) 

CX:iSD('n/6-311+G(2dfl 

//QCISDCn/6-31 1  +G(2df) 

//QCISDCn/6-31 1+G(2df) 

//QClSDCn/6-31 1  +G(2df) 

Rs(Na-Ai)s3.182A 

Ra(Na-Ai)=3.007  A 

R«(Na-Ai)=2.739  A 

EmP4=*403.79651 

EPMP4=-403.79027 

EpMP4=-403.76274 

EQClSDm=-403.80002 

EQCISDm=-403.79176 

EQCISDrn=-403.76424 

<s2>=2.027 

<S2>=2.018 

TeMP4=0.0  kcal/mol 

TePMP4=3.9  kcal/mol 

TePMP4=21.2  kcal/mol 

TeQCISDm=0  0  kcal/mol 

TeQCISDm=5.2  kcal/mol 

TeQCISDm=22.4  kcal/mol 

Table  X.  Calculated  Molecular  Properties  of  the  Lowest  MgAI  states. 


Table  XI.  Calculated  molecular  properties  of  the  lowest  LiAl  states. 


Table  XII.  Calculated  molecular  properties  of  the  lowest  MgSi  states. 


MgSi  (3nr) 

1 0^20^1  It 

CASSCF-MrtCISD 

CASSCF-MRCISO 

Re(Mg-Si)=2.S56  A 

R«(Mg-Si)=2.768  A 

(0e=317  cm’’ 

©6*271  cm*' 

EcASSCF-MRCISD='480.62O1  1 

EcaSSCF-MRCISD=-488.591 48 

T8casscf-mrcisd=0.0  kcal/mol 

TeCASSCF-MRCISD=18.0  kcal/mol 

CASSCF-MRCISDfQ) 

CASSCF-MRCISDfQI 

R«(Mg-Si)=2.5S9  A 

R6(Mg-Si)*2.876  A 

<b«=314  cm*'' 

©6=264  cm*1 

EcASSCF-MRCISDfQ<=-488.62663 

ECASSCF-MRCISD/Q)=-488.60041 

TeCASSCF-MRCISDfQi«0.0  kcal/mol 

T6CASSCF-MRCISDrQ1»16.5  kcal/md 

Table  Xlli.  Calculated  molecular  properties  of  the  lowest  BeP  states. 


Table  XIV.  Calculated  molecular  properties  of  the  lowest  BP  states. 


BP  {3ni) 

BP{ll+) 

1<j22a^lR^3ot 

1 0^20^1 

CASSCF-MrtCISD 

CASSCF-MRCISD 

Ri,(B-P)=1.755  A 

R«(B-P)a1.684  A 

(aas920  cm‘^ 

<Be=t038  cm’^ 

EcaSSCP-MRCISD=-365.53221 

EcASSCF.MRCISD*-36S.52083 

TaCASSCFMRCISD=0.0  kcaVmol 

TeCASSCF-MRCISD=7.1  kcal/mol 

CASSCF-MRCISO(Q) 

CASSCF-MRCISCHQi 

Ra(B-P)=1.765  A 

Rb(B-P)»1.687  a 

(i)a=897  cm"^ 

o)e=1026  cm'^ 

EcaSSCF-MRCISD/O>='365.55024 

EcASSCF-MRCISDIQ\=-365.53602 

TeCASSCF-MRCISD/Q)=0.0  kcal/mol 

TaCASSCF-MRCISDIQ)=8.9  kcal/mol 

-> 


Table  XV.  Calculated  molecular  properties  of  the  lowest  SiP  states. 


HF,  11+  1 

HMg.  2] 

r+ 

HAI,  11 

HSi.  20 

r 

HP.  31- 

HS.  20i 

HCI,  11 

U2,  iSa 

LiBe,  22 

'+ 

Table  XVI.  Calculated  and  experimental  ground 
of  diatomic  molecules  contains  first  and  second 


Molecule  D  (eV)  calculated 


HoJIa-"  4.47807ia 


2.415c 


2.047c 


3.6470 


3.443c.  3.5770 


HN.  35;-  3.339c.  3.4996 


4.380C,  4.5180 


5.915c  .  5.9840 


1.27 


3.1749 


3.1619 


3.1359 


3.7219 


4.6539 


1.02,  1.1 40C 


0.26^’,  0.29',  0.251 


LiB,  311  1.18^.1.11 


LiC,  2.58 


UN,  3X*  1.49",  1.61 


LiO,  20'  3.48"  .  3.78',  3.30 


LiF,  il-*-  5.81"  5.980C  ,  6.06 


LiNa,  il-'-  0.85 


LiMg,  22;+  0.18",  0.20' 


LiAl,  ■'I-'-  1.01 


LiSi,  1.83 


UP,  31-  1 .66 


US.  201  3.08",  3.30 


LiCI,  ■'1+  4.76"  4.89'.  4.86 


Be2,  ■’la'"  0.09'.  0.08-0.10° 


BeB,  20  1 .579 


BeC,  3S-  2.39^ 


BeN.  1.34s 


BeO,  4.69 


BeF.  2S+  5.94',  5.71 


BeNa,  2X+  0.14' 


BeMg.  11+  0.008‘ 


BeAl,  20  0.40' 


BeSi.  3X-  1.28" 


BeP,  1.06^ 


BeS,  11+  1 2.29',  3.23 


state  dissociation  energies  (D) 
row  atoms 


D  (eV),  experimental 


4.478077 


2.429 


<3.47 


4.392 


5.869 


1.34 


<3.06° 


3.1709 


3.06°  ,3.1619 


3.1229 


3.55°  .  3.7649 


4.434°  ,  4.6539 


1.046 


1.54m 


2.53" 


4.60 


5.85°,  6.26 


BeCI,  21+ 


Bo.  31 


BC. 


BN.  3n 


BO.  21+ 


BF.  1S+ 


BNa,  3n 


BMg. 


BAI,  31- 


BSi,  41- 


BP.  3n 


BCI, 


Co,  lla^ 


CN.  2X+ 


CO.  11+ 


CF.  2n 


CNa,  41- 


CMg. 


CM.  4^- 


CSi,  30; 


CP,  2S+ 


CS.  11+ 


NO.  2ni 


NF.  31- 


NNa, 


NMg. 


NAI,  30 


NSi,  21+ 


NP.  11+ 


NS.  20 


02.31 


OF.  20 


ONa.  20 


OMg,  11+ 


OM.  2S+ 


OSi.  11+ 


OP.  20 


OS.  31- 


OCI,  20 


II 


3.87'.  3.84 


2.84W 


4.21* 


4.56y 


8.43^.  8.32 


7.74 


0.76^ 


0.47q 


1.78^ 


3.15^ 


3.13^ 


5.71 


5.33f.  5.49 


6.26^  6.40 


7.623C  . 


11.231C 


5.71 


1.97q 
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